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Abstract. In this article the quasi-Keplerian parameterisation for the case that spins and orbital angular 
momentum in a compact binary system are aligned or anti-aligned with the orbital angular momentum vector is 
extended to 3PN point-mass, next-to-next-to-leading order spin-orbit, next-to-next-to-leading order spin(l)-spin(2), 
and next-to-leading order spin-squared dynamics in the conservative regime. In a further step, we use the expressions 
for the radiative multipole moments with spin to leading order linear and quadratic in both spins to compute radiation 
losses of the orbital binding energy and angular momentum. Orbital averaged expressions for the decay of energy 
and eccentricity are provided. An expression for the last stable circular orbit is given in terms of the angular velocity 
type variable x. 
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1. Introduction 

Recent progresses in the post-Newtonian treatment of compact binary systems with spinning components call for 
an extension of the known parametric solutions to the binary dynamics to include latest spin interaction terms as 
well as radiative dynamics. As computer resources are recently unable to generate thousands of orbits, which are 
interesting to ELISA (which may be, optimistically, able in several decades to see gravitational waves of orbital 
periods of a few hours), it is desirable to have accurate and efficient gravitational wave templates for an analysis of 
the detector data. In a prequel paper [1] reasons can be found why to regard compact binary systems with aligned 
spins with regard to their implications to gravitational wave data analysis. This article will aim to incorporate 
higher-order terms to the in the orbital and spin dynamics to the quasi-Keplerian parameterisation for compact 
binaries with aligned spins. These may become interesting when the binaries come to the final stage of their live 
before merger. Let us, for convenience, state some milestones in the recent literature. For literature on point- 
mass Hamiltonians through 3PN and leading-order spin-squared and spin-orbit Hamiltonians, we refer the reader 
to the introduction of our previous paper [1] and start from there. Note that there were two recent publications 
concerning the 4PN conservative point-mass dynamics (see [2] for the 4PN Lagrangean and [3] for the Hamiltonian 
in the center-of-mass frame, both as preliminary results up to order G^). 

The next-to-leading-order (NLO) spin-orbit (SO) contributions have been derived in [4] and later in [5]. Next- 
to-next-to-leading order (NNLO) spin-orbit and spin(l)-spin(2) (S(1)S(2)) Hamiltonians are recently derived in 
[6, 7] and the corresponding Lagrangian potential via the effective field theory formalism in [8] and will also be 
used in our calculation. The leading-order spin(l)-spin(2) Hamiltonians are available in [9] and extended to next- 
to-leading-order in [10]. Spin-squared dynamics (S(l)^, S(2)^) depend on the model - or more precisely: on the 
equation of state - of the matter of the constituents of the binary. Depending on its rotational velocity and its 
stiffness, the included body will deform and self-induce a spin quadrupole moment that will start interacting with 
the binary orbit and re-couple to the gravitational field. The proportionality factor Cq will represent this issue. 
It varies from one of black holes to four for neutron stars (such that it is related to the constant Aq of object a 
in Ref. [1] via Cqa = — 2Aa) and thus characterises how a body resists rotational deformation. References [11, 12] 
provide the next-to-leading-order interactions of this type. 

Anyway, as the spins do not precess, the quasi-Keplerian parameterisation [13, 14] can be employed to obtain 
a parametric solution to the dynamics. When the spins are not aligned, their orientation is not constant. There 
exist several publications about precessing spins, see e.g. [15] for the case of "simple precession" (which means 
circular orbits together with the fact that the angle between the total angular momentum J = L -|- Si -|- S2 and L 
is conserved), [16, 17, 18] for the case that only one body is spinning or the masses are equal, and [19] for circular 
binaries with unequal masses. 

As gravitational waves carry away energy and angular momentum, the semimajor axis and the orbital 
eccentricity will suffer a slow decay in the validity regime of the post-Newtonian approximation. The radiative 
losses of compact binaries have been extensively discussed in the literature. Reference [20] gives a general expression 
for the losses due to gravitational waves in terms of the mass and current-type multipole moments of the binary. 
This has been elaborated in general in [21] and applied to nonspinning compact binaries in [22] through 2PN point 
particle dynamics, and, recently, in [23] to 3PN point particle dynamics. In Section 4 we give expressions for the 
decay of the orbital energy and the radial eccentricity, deduced from the spin dependent multipole moments given 
in [20] and [24]. 

Since we used different approximation levels for the last stable circular orbit calculations, the quasi-Keplerian 
parameterisation, and the radiative dynamics they will be applied exclusively in the corresponding sections. This is 
due to the fact that the radiative and dissipative effects are known for the spin up to different approximation levels. 
The next-to-next-to-leading order spin-orbit Hamiltonian completes the knowledge of the dynamics of binary black 
holes up to and including 3.5PN for maximally rotating objects. For general compact objects like neutron stars 
the leading order spin(l)^ Hamiltonians are still missing. Since the NNLO spin(l)-spin(2) Hamiltonian is at 4PN if 
both objects are rapidly rotating, the full post-Newtonian approximate dynamics up to and including 4PN requires 
further efforts. We state that for general compact objects the leading order spin(l)^ Hamiltonian is only known for 
black holes and has to be generalised to include neutron stars as well. Furthermore, the NNLO spin(l)^ Hamiltonian 
- which is also at 4PN if the object is maximally rotating and maybe even stronger than NNLO spin(l)-spin(2) - 
is completely unknown. 

The structure of the paper is as follows. First of all the determination of the last stable circular orbit and a 
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resummed binding energy will be discussed in Section 2. Afterwards in Section 3 a quick summary of the quasi- 
Keplerian parameterisation for eccentric orbits and the appropriate orbital elements will be given. In Section 4 the 
dissipative dynamics and energy and angular momentum loss will be discussed. Finally, in Section 5 the conclusions 
and future applications will be provided. In the appendices there is a table with symbols and a short subsection of 
the rescaling of several quantities to simplify equations and discussions. 



2. Last Stable Circular Orbit 



To determine the last stable circular orbit our starting point is the binding energy for circular orbits. Either one 
can extract it from a Lagrangian potential constructed from a given metric or one can get the binding energy from 
a Hamiltonian in the center-of-mass frame for circular motion. This means in the center-of-mass frame Hamiltonian 
the Pr = (ni2 • p) component of the linear momentum must be set to zero and = L^/r^ holds. For a system with 
spins we further set the spins to a configuration in which they are aligned with the orbital angular momentum. For 
an energy E{r,L) the circular orbit is given by the perturbative solution of 

^(r,i) = 0, (1) 

for L. Then E{r, L{r)) is the binding energy at the circular orbit r. Since the radial coordinate depends on a certain 
gauge we have to transform it into a gauge invariant quantity (invariant up to a change in the spin supplementary 
condition) given by 

\ 2/3 

, (2) 

which is related to the orbital angular velocity, see Appendix A and Appendix B. For the Schwarzschild spacetime 
the binding energy for a circular orbit is given by 

1 — 2a; 

£^Schw(a;) = ^ -I. (3) 

For a test-spin in a stationary Kerr spacetime we can also write down this expression for orbital momentum aligned 
Kerr spin and test-spin to linear order in test-spin (see [25, 26, 27] for the appropriate potentials^), namely 

2 

£^KGrr,Ts(a;) = — 1 ^ , ^ 

V — a + x^-^/'^y a + x^^/^ — 3v— a + x^'^l'^ 
3.-3/2 

-(- 



V-a-|-a;-3/2-v/a-|-a;-3/2 - 3v^-a + .T-57^ 
( -a + \J-a + a;-3/2 ] 2;3/2 

+ S \ ' . (4) 

V-a + a;-3/2Va + x-'^l'^ - 7,^ -a + x-'^l'^ {ax^l"^ - l) 

In the post-Newtonian case one can write down this binding energy also, but we will not provide it here since 
the expression is very lengthy and gives no further deep insights into the calculation (see e.g. [28] for very recent 
results). We briefly note that it is a polynomial in yfx which corresponds to a post- Newtonian expansion. From 
the binding energy for circular orbits E{x') we constructed a quantity e{x) via 

T 2 



e{x) 



(rjEix) + ly -1 ^ 



-1, (5) 



2r/ 

(see e.g. [29]). We refer to e{x) as "modified binding energy" in contrast to "binding energy" in case of E{x). The 
modified binding energy e{x) has for a test-mass (77 — 0) moving in a Schwarzschild spacetime a polynomial structure 
in X in the numerator and denominator. For the binding energy of a test-spin moving in the equatorial plane of a 
stationary Kerr black hole this is not true anymore. The relation between the xi, X2 spin magnitudes and a, S is 
given by xi = a and X2 = S/j]. Notice that the l/r] terms are not singular, because 1 — y^l — Arj = Arj/(\ + ^l — Arf), 
which renders these terms well-defined at 77 = 0. This issue appears due to the fact that S" is a test-spin and so may 



X In the given literature tliere are a few typos which were corrected in the appendix of [28]. 
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not vanish in the Umit 77 — > 0. Here one has to approximate in Kerr spin and test-spin to get a rational structure. 
The mentioned modified binding energies are given by 

l-Ax 



eschw(a;) 

eKerr,Ts(a;) 



= — X 



l-3x 



(6) 



-^-a + a;-W/2 + x^{a^ - 3a</-a + x-V^ + 4^-a + x-^/^^) 
(-1 + oa;3/2) ^-1 + a:3/2 + 3^-o + x-3/2)^ 

(-a + ^-a + a;-3/2) (-1 + 2v'-a + x-3/2a;3/2) 



25- 



(-1 + aa;3/2)' (-1 + [-a + 3v^-a + a:-3/2] 3,3/2^ 



In the approximation in Kerr spin and test-spin the modified binding energy reads 

^approx.Kerr,Ts('^) ~ 



1 - 4x + |aa;3/2 _ ^2^2 ^ |^^5/2 _ 1^^2^3 ^ 4^2^4 

l-3x + 4ax5/2 _ o2.x3 - |a2a;4+25'.x-V2 - 2aSx3 + 2S'xV2 - jaSx^ 



(7) 



, (8) 



where one can still identify the terms coming from the test-mass motion in a Schwarzschild spacetime. (Notice 
that the test-spin parts in the denominator were implemented by using the geometric series at first order, i.e. 
1/(1 .t) w 1 + .X, to implement the correct pole structure coming from the test-spin.) This rational function in 
y/x can be taken as a starting point interpolate between test-spin modified binding energy and post-Newtonian 
approximate gravitating mass orbiting another gravitating mass. 



2.1. Construction of Binding Energy 

After having an initial guess for the rational modified binding energy we tuned all parts of the numerator by 
rj-dependent coefficients and matched the Taylor expansion in ^/x with the post-Newtonian approximated e{x) 
obtained from the ADM-Hamiltonians (3PN point-mass, ncxt-to-ncxt-to-leading order spin-orbit, ncxt-to-ncxt-to- 
leading order spin(l)-spin(2), next-to-leading order spin-squared). Since we only had linear corrections in the 
test-spin to start with, we tuned the denominator by adding linear-in-test-spin contributions with the opposite 
sign to get the correct singularity structure of the modified binding energy e{x). These considerations lead to the 
expression 



^tunedKerr,TS (^) " ^ 



+ a 



S 



10 r. ^/43r/ 



.3/2 



91t] 
IF 



47?? 

u ~ 

4ry2 
9 
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"9" 
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81 



55 43 n 2 , r— 

18 + 18 ^^+3^(^ 
2057r2r7 36797? \ 3 



36 



+ a 



Cqi 

95r?2 

^6^ 
91tj 
6" 

+ ( F(r/) + aS 
4 



_4jf 
3 

197r/ 



+ 



96 



72 J 
7\ _ 10772 
6/ 9^ 
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577 
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18 9 
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Q 1 

+ 3^/^^+^(l-^/^^) 



+ Za^frirj)]x' 



1 - 3a; + (4a + 2S)x'^/^ - (a^ + 2aS)x^ + 2Sx'^'^ - , («' + ('S)x'^ 



(9) 
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Figure 1. Last stable circular orbit for S = 0.1 plotted for different symmetric mass ratios rj and Kerr spins a 
(—1.0—0 from blue to black, 0—0.8 from black to red, 0.8-1.0 from red to violet). Notice that the approximation 
is not valid anymore above the dashed line. Also notice that to large Kerr spins a are invalid because the last 
stable circular orbit will be of the order of magnitude of the Schwarzschild radius which violates the post-Newtonian 
approximation (wide separation). 



There are no quadratic terms in the test-spin, so there will be no Cq2 in the expression. However, S* is a test-spin 
and the denominator of e(x) is only valid in the test-spin limit, hence the results given here are only valid around 
77 « 0. The unknown function F{r]) has to be fixed by the 4PN point-mass Hamiltonian later (see [2] and [3]) 
and the unknown function f7{r]) by the next-to-next-to- leading order Sf Hamiltonian ifriri)^^^^! is required to be 
consistent with the Kerr-limit). 

2.2. Last Stable Circular Orbit 

The equation detuncdKcrr.Ts(a;)/da: = is numerically solved and one can obtain the solution .tlso for certain mass 
ratios r], spins and quadrupole constants, see figure 1. 

To link to recent literature, we wish to mention reference [29] where the last stable circular orbit through third 
post-Newtonian order for point masses has been computed. To generalise to certain configurations including spin, 
the binding energy of "last stable spherical orbits" has been derived, for example, in [30] in the effective one-body 
approach for non-aligned spins. Compact binaries with spin under next-to-leading order spin-orbit coupling evolving 
in circular orbits were studied in [31]. We also wish to reference [32] where corrections to the last stable circular orbit 
of a Schwarzschild black hole due to the gravitational self-force have been derived, and also [33] where the authors 
calculated gravitational self-force corrections to strongly bound eccentric orbits in a Schwarzschild spacetime. 

3. Eccentric Orbits: The Quasi-Keplerian Parameterisation 



The next-to-ncxt-to-leading order spin-orbit, the next-to- and next-to-next-to-lcading-order spin(l)-spin(2), and 
finally the next-to-lcading order spin(l)-spin(l) Hamiltonian in reduced form (in the center of mass) are listed at 
the beginning of this section. The point particle and leading-order spin Hamiltonians can be found in [1] . We define 
the sums and differences of the two canonical spin tensors as 



3.1. Overview 




Si + S2(i)(j) , 



(10) 

(11) 
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labeled with a "hat". Those satisfying the covariant spin supplementary condition will be labeled with a "tilde", 
namely 



(12) 

(13) 



and the latter will become important for the multipole moment expressions taken from the literature. The reduced 
Hamiltonians read 



H. 



NNLO _ I 1 

SO - 



- 8^) (P • ni2)' {n\2P'^ij) + (nl^p'^j) (^(45r/ + 256)(p • n.^f 



- ^(39r, + 314)(p2) j + ^(737, - 206) (nl^p'S,,) (p^) 



-15 



r?2(277-l)(p.ni2)Mnl2P^S 



3 

32 



H 



NLO 



1 



Si 



V i&rf - Uv + 4) (p . n,,r {n\,p^±,j) (p^) + ^/^^ {n\,p^\,) (^lME_i^ 
r?(9r? - 4)(p • ni2)^(p^) + ^7?(22,7 - 9)(p^)') - {^v' -877 + 2) (ni^P^E,,) (p^)' 
^ (-2772 + 337? + 42) (nl^p^tij) + ^ v'r^(7y + 1) (nl^p'Aij) 



S(1)S(2) - ^4 ] ^3 



-y??^(p-ni2)^(ni2n^25'i(j)(q)S'2(fc)(g)) + ('S'i(j)(fe)S'2 0)(fc)) y (3(p • riia)^ + (p^)) 



- ^'?'("i2"^2'5lO)(g)52(fc)(g))(p2) + \{p'p''Si^j)(g)S2ik)ig)) 

+ (P • ni2) ((2t? + 1)?? + Vl - 4w) (n^Si (i)(q)P^ S2 o)(<,)) 

+ (»?(2?? + 1) - - (»^i2'S'2(i)(fc)P^'S'2(j)(fe))^ 

- ^(^ + 2)??(ni2F''S'i(i)0)) (7ii2p''S'2(i)0)) 



^("i2"l2'5lO)(g)'S'2(fc)(g)) " 3r? (S*! 0)(fc) S'2 0)(fe) 



) . 



H. 



NNLO 



SI 



s(i)S(2) - c6 1 r-i 



+ 



77 ^7l'{2^12'^l (j)(g)'S'2(fe)(g)) (303(p • ni2)^ + 77(p2)) + — 7?^(p'p''S'l(j)(,).S2(fc)(g)) 

\il (l2(47? + 3)(ni2P^-^i(i)(,)) (nl2P^'^2W0)) -r?(^i0)w52 0)(fe)) (52(p • ni2)' + 23(p2))) 



+ ("l2^1(i)(9)P''^2 0)(g)) - 4W(?? + 3)(P • ni2) - ^7/(447? + 9)(p • ni2) 

+ {r)^i2S2(i)(k)p'S2U)(k)) 0Vl-4w(r; + 3)(p • ni2) - ^7/(447? + 9)(p • ni2)^ 

(nl2^i(i)(,)P'^2 0)(,)) (^r;(p • 1112) (207;2(p • ni2)' + + 4,? - 3) (p^) 



+ 



+ 



16 



Vl - 4777?(277 - l)(p • ni2)(p')j + (ni2^2(i)(fc):p'52 0)(fc)) (^7?(p • ni2)(207?2(p • 1112) 



(14) 



(15) 
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NLO 

S(l)2 



ttNLO 
-•^3(2)2 



(677^ + 477 - 3) (p2) j + -^/T^(l - 2r7)77(p • n,2){p^] 
^772(pV5i(,)(,)^2(fc)(g)) ((4r; + l)(p2) - 6r/(p • ma)' 



16 
1 



r,f 30r,2(p . ni2)'(5i(,)(,)^2a)(fe)) + (p") (277(87/ - 2)(p2) (^^ (^-)(,)^2 



3 (477^ + 177] - 6) (77.*2P^S'i(,)(j)^ ("l2P^'S'2(j)0) 

677(p • 1112)^(577(7112^ S*! ) { 71^2^ 5*2 WO) ) - 2(7/ - l)(p^) { Si (j)(fe)52 (j)(fe) ) ) 

3 



SI 



(<2 

772(ni2nt2^iO)(,)^2(fe)(g)) (3577(p • ni2)' + 10(77 - l)(p • ni2)'(p') + (877 - 2)(p2)') 

^77(1077 + 7) (s'i(j)(fc)S'2(j)(fe)) - ^?7(7977 + 105)(7^{2"^2'S'lO■)(g)'S'2(fc)(g) 

/ „ ^ \ / 15 15 \ 

(p • nl2)^(7^j27^^2-5'l(i)(9)S'l(fc)(9)j ( yC'Qi77(277 - 1) - y Cqi^I - 47/77 j 



(16) 



+ (^1 (MkAuKk)) (^(P ■ ni2)' (-^'7(Cqi(477 + 2) + t; - 4) - ^(Cqi - 2)^1^77^ 
+ (p') Qyr^((CQi - 2)77 + Cqi) + ^(Cqi(6 - 27/(277 + 3)) + 3(7/ - 4)77)^^ ^ 
+ (p • ni2) (7ii2^i 0)(,)/^i (fc)(,)) Q^(-2Cqi(7/ - 2) + 7/ - 3) + ^(4Cqi - 3)^1^7/ 
+ ("i2«^2^i UKiA (km) (P') (^(3Cqi(77 + 3) (27/ - 1) - 3(77 - 4)77) - ^yr34^(CQi(7/ + 3) - 47/)) 
+ (pV^I (,)(,) ^l(fe)(g)) (-|r/(3CQi +77- 3) - ^(Cqi - 1)^1^7/ 

(^i2'^i2^iO)(,)^i (fc)(,)) (J(3Cqi(77 + 1) - 167/ + 6) + i/r34^(CQi(97/ + 3) - 47/ + 6)) 
+ (^1 UKkAum) {li-CQiiv + 1) + 37/ - 1) + i/r^(-CQi(37/ + 1) + 7/ - 1)) 



H^(LO (1^2), (18) 

where Equation (18) follows from the fact that ni2 always appears in a quadratic form and the sign has no influence. 
The quasi-Keplerian parameterisation is the basis of the calculation of the radiation losses. Having polar coordinates 
(r, 0) for the plane of motion characterised by L and the spin as = Xa^/L for the ath object at initial instant of 
time, we list the elements of the parameterisation schematically without stating technical details of the computation. 
These can be found in, e.g., [13, 14] and [1]. 

r = Or (l — Br cos £) , (19) 

Af{t - to) = 5 - et sin f + F^_s {v - £) + sin v + sin[2v) + Fg^, sin(37;) , (20) 
2ti , 

- (/)o) = w + G2v sin(2w) + G^v sin(3w) + G^^ sin(47;) + G^^ sin(57;) , (21) 



:= 2 arctan 









- tan 







(22) 



It is needed to re-express the luminosity and angular momentum loss which are given in general terms of v^, and 
r to get a more or less explicit time dependency (note that the Kepler equation (20) cannot be inverted by t — 
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exactly). 

3.2. Semimajor axis and eccentricities 

Let us define en as the "Newtonian" value of the orbital eccentricity, 
CN y/l-2\E\L^. 

Then, the orbital elements read as follows. 



(23) 



2\E\ c2 I 4 
+ 



+ 



^ I (xi - X2) + (i - 1) (xi + X2) I 

^ I iCqi (2,7 - x/r^ - l)x? + \cQ2i2r^ + x/T^ - l)xl - VX2X1 1 
(XI - X2) (s - ^) ^ + (XI + X2) 



1. 

C4 1 L3 



8L 



(Xi - X2)(5?7 - 8)^1-4??+ (-6r?2 + 19r? - 8) (xi + X2) 



{ Y [i L^Q^ (4 + - + (t - 2) ^^^) - 2 - 9r? + 3) + (6 - 6^)yr^) 



+ 1^(cq, ((,-1)^ + 0:34,) 



+ 



+ 



4^*^' {\ ^'^"^^ + ^^'^ ~ + ~ ?) ^^^) - 2 (ry' - 9ry + 3) + (6 - 6r?) 



xi/ 2 

2 

m 

L2 



3r]^/l^^ 1 



(77 - 3)77 



?7(t? + 32) _ (77-4)77!^! 

L4 L2 



X1X2 



1 I 1 

7s 



+ 1(77^-2772- 377 + 



+ J I ^ (^^1^^7(1377 - 3) (xi - X2) - (1 V - 2I77 + 5) (xi + X2)) 

+ ^ (1277' - 12977^ + 3I877 - 130) (xi + X2) - (2977^ - 1497, + 130) (xi - X2) 

+ ( x/r^ (477^ - 6977 + 97) (xi - X2) + (-v' + ^ - ^ + 9?) (xi + X2)) | 

77|E| 



^1 



77_W 
2L2 



+ 7t X1X2 <; - ('^ - 1) - (2 V + 191?? - 373) + ^ (1377^ + I2577 - 581) ^ , 



(24) 
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1 L2 



+ L2(17 _ 7r^)\E\^ + ^ (2{n - 2)(xi + X2) - 4./T^{xi - X2)) + 4(7? - 1)\E\ | 

+ (-16772 + 4777 - 112) \Ef + 1^(17 - llT?) + ^^^^\/2eg,(277 - 5) + 2 (St?^ + r? + 2) \E\'^ 



+ Ss 



4^2 



eg, (772-87? + 6)|E|3/2(xi+X2) 



+ 16v/2^(77 - 3)|i?|5/2(xi - X2) 

+ \E\ (^^ ((-4772 + 3977 - 32) (xi + X2) + yr^(977 - 32)(xi - X2)) 

" (^/-3)(xi-X2)\/(2-87?)|i;| 
2 1597? 



+ 



L2 



( (l6??' - ^ + 62^ (xi + X2) + ^0^(124 - 597?)(xi - X2)) 



+ 



§ I ("^^1 ('^'^^^N (-6??' + 33r? - 14) + eg, (-4r?2 + 29r? - 8)) 

+ 4v/2xiX2e^(7? - 15)7? + V2x^ (CQ2e^ (-67?2 + 33r, - 14) + (-4r;2 + 29r? - 
+ ^ (^Xl (-y<^Qi (V - 777 + 3) + (y<^Qi(^ - 3) + 26?7 - le) - 9772 + 587? - 16^ 

+ xi (-y C'Q2 (V - 777 + 3) + (-y Cq2(7? - 3) - 267? + 16^ - 97?^ + 58?? - 16^ 

+ 2xiX2r?(107? - 73) j 

+ \E\ (jj xl (Cqi (67?2 - 4l7? + 18) + yr^(CQi(18 - 57?) - 24(77 - 1)) + 8 (77^ - 97? + 3 
+ xi {Cq2 (67?2 - 4l7? + 18) + yi - 477(Cq2(577 - 18) + 24(r? - 1)) + 8 (7?^ - 97? + 3)) 

+ 4xiX2??(»7 + 32) 

+ 7I [ X?(Cqi(57; - 14) + 137? - 8)^\E\-47i\E\+ xi(CQ2(14 - 57,) - 137? + 8) - 477|i^| 
+ (x?V2-8r?(CQi(28 - IO77) - 267? + 16) + xi ^2 - 8r?(CQ2(107? - 28) + 267? - 16)) | 



+ 



c6 1 i4 



\El 

i2 



(I6772 + (577^ - 281) 7? + 134) 



|^|3/2 

24^22.3 



(14407?^ + (I237r2 - 13952) 7? + IOO8O) 



I67772 /'3127 397r^ 



+ 



r,-^\+L' (-SOr,' + - 250, + 66o) li?!* 
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- 3V2L (467?2 - 193?7 + 265) \Ef/^ + (324077^ + {Wdn^ - 19964) 77 + 16380) 

+ (20?7^ + 4277^ - 13577 + 46) \Ef\ 



+ 

c 



i^sj ^ Qyr^ (171??' - I2II77 + 1384) (xi - X2) 
+ (-34773 + ^- ^+692) to+x.)) 

+ ^-^^ (^2 (4l7/-^ - 665772 + 18447? - 966) (xi + X2) - 2^2 - 87; (50r/2 - 3947? + 483) (xi - X2) 
+ \E\ (^^ [(47?3 - 1587?2 + 6507? - 388) (xi + X2) - 4^1 - 47? (47?^ - 697? + 97) (xi - X2) 



^ {7v- - 997/ + 140) (xi - X2) V(2 - 877)1^1 



L4 



+ ^ (- V + 6977^ - 23977 + 140) \Ef/\xi + X2) 



|3 



? ( ('^^''' " + ^^^"^ ~ ^^^) Vl-4?? (2437?2 - 77O77 + 996) (xi - X2)) 

+ \E\y^ (^4^2 - 87? (277/2 - 1337/ + 177) (xi - X2) 

+ 2V2 (-277/3 + 2777/2 - 6507? + 354) (xi + X2) 



1 .2 j 2V2X1X2V (15??2 - 44277 + 939) lEl^^ X1X2V (407/2 - 7517/ + 1364) 

n/2xiX2?? (-517/2 - 13067/ + 38 49) \E\^^'^ IOX1X2V (H^^ + 1367/ - 354) \E\'^ 
6V2xiX2r/ (7r/2 + 83r/ - 315) 4xiX2r/ (137/2 + 125r/ - 581) \E\ \ 



+ 



2 _ 2 



i6 



+ L 



+ 



+ + L'ifj - 15)\E\^j + (s^T^ixi - X2) + (8 - 47/)(xi + X2)) 

^ (4(7/ - 2)(xi + X2) - 8vT^(xi - X2)) I 
§ I ^ (x? (Cqi(3 - 67?) + 3Cqi v/r^) + 12xiX2?7 + xi (<^Q2(3 - 677) - ^0(^2^1^)) 
\E? (xl (Cqi(87/ - 4) - 4CQiyr^) - 16X1X2?? + xl {Cq2{8v - 4) + 4Cq2^/T^) 



+ A + 15. - 80^ 1^1^ + W f-^ - 297/ + 5lU f ¥ + 447/ - 8 ) \Ef 



2 ' 8 ' y V2 

+ J {(Xl + X2) (-|^ iv' - 71v + 64) + ^ (7/2 - 367/ + 17) + L(80 - 3l7?)|£;|3) 
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+ (Xi - X2) (^1^^1^(1177 - 64) - L^/l^{rj - 80)\E\^ + ^(34 - ISr?)^!^ 
^si 2f\E\(„ f^W 1261r? „„\ , — 205r7\ 291?] 



^ ^+72 Cqi(-30,7 +221,7- — 



+ v/l-4r/ ( Cqi ( 50?7 - ^) + lOS?? - - ^ + 215r? - 56 



+ \E\^ (^Cqi (iv^ + - 28^ + (^Cqi (^-^ - 28^ + 9r? - 48^ - ^ + 105,? - 48 



+ 1^-5^+72 



Cq2 (^-30r/2 + 22177 " ^) + ^1 " 4r? (^Cq2 - 50??^ - IO377 + 56 



+ 215?? - 56 j + |£;|3 (^Cq2 (^7772 + ^ - 28^ + y^I^ (^Cq2 (^^ + 28^ - 977 + 48^ 



697/^ 

^ + 105. - 48^ + V, V, ^ ^f857, + 1484) - 



2 + 1057? - 48j J + X1X2 (^^(85r? + 1484) - ^(9r? + 290) - 7?(197? + 292)|E|3 j | 
+ ^ I (-707?^ + 34407/^ + (I3257r2 - 65436) 7? + 27776) - ^L^ (77^ + 137?^ - 807? + 248) 

+ (-107?^ - 26 897?=^ + (26860 - 58l7r^) 77 - 2456) 
+ ^ (234r/^ + 45 367?^ + (2764 + 37r2) 7/ - 16032) \Ef^ 
, '^s f, J \E\ 



(Xi - X2) (- ^ V^T^ (16177^ - 79 687? + 10336) 



\E 



+ ^^v/l-47?(l2977?2 - 279887? + 13568) 

+ ^L^1-4tj (3I772 - 2147; + 2152) 1^1" + (-237772 + 404?? + 3664) 

+ (1^ (_ H|! _ + 2.78, - 1292) + ^ ^ 7^ ^ 20^ _ 12041, ^ 

+ L (-2,3 + _ + 538) + H! (^ - - 737, + 916 

+ ^ I {-^L%r]{63r] + 220) + 5728)|£|3 - 8i^(5??(5?? + 107) + 5247)|i;|^ 
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+ 5L^{T]{335r] + 5306) - 8552)|^| - 7567?^ - UdWrj + 61766) I , (26) 



e2 = + 1 {{rj - 3)\E\' - 2{rj - 6)\E\} 



+ S { ((X. - X.) {sLVT^m^ - ^v^lZffl) ^ + (,(s - 4,)|^|^ + li^) 
^ |^|(C,.(2-4,H2C,.^) ^ (^^^^^^ _ _ J 

+ 1 (- V + 55,7 - 80) \E\' + (34-22^)1-^1 + (^2 ^ ^ ^ 26) | 
+ J I (Xi + X2) (-8r?' + 78r? - 64) + i {6r,^ - 49r? + 80) \E\^+^-^ (IOt?^ - 70r/ + 4) 

+ (Xi - X2) (^^1^(97? - 32) + L(80 - 19v)Vl^\Ef + ^(1 - Arjf/"^^ | 
+ § r M (Cqi (I2r?2 - 82r? + 36) + ^/T^{CQ^{36 - IQrf) - 48(r? - 1)) + 16 {rf - 9?? + 3)) 



+ ^ (Cqi (-187?2 + 85r? - 37) + ^J I - Ar]{C Qi{llr) - 37) + 52,7 - 8) - 207?^ + 68?? - 8) 

+ \E\'' (Cqi (-4772 + 34r/ - 14) + ^1 - 4r/(CQi(67? - 14) + 12(77 " 4)) - 2 (r;^ - 54r/ + 24)) j 
+ ,,,, (?2(II±M _ ^ _ ,o)„^|s) 

+ X2 (Cq2 (I2r?2 - 82r? + 36) + ^1 - 4r?(CQ2(107? - 36) + 48(r? - 1)) + 16 [rf - 9?? + 3)) 

+ ^ (^Q2 (-187?" + 8577 - 37) + ^1 - 477(Cq2(37 - II77) - 5277 + 8) - 2077^ + 6877 - 8^ 

+ \E\^ (Cq2 {-^V^ + 3477 - 14) + V1 - 477(Cq2(14 - 677) - 12(7? - 4)) - 2 (77^ - 5477 + 24)) j | 

+ i m (32,^ + 2 (5.^ - , + 268) + ^ (-40,^ + (?f^ - 5^) , - « 



+ 3897? - 372^ \Ef + ^-97?^ - ^ (l72 + 37r^) 77 - 32^ \Ef^ 

(5s f / 8V1 - 477 (4772 - 6977 + 97) \E\ ^T^^ {iQif - 5977? + 352) \E\'^ 

+ ^LvT3477 (3577^ - 58177 + 1076) li.^ - + ^^^^ " ^^^) "^"^^ 



2 V - V - -./II 2L 
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, / 4 (27]^ - 79r?2 + 32577 - 194) \E\ (-20r]^ + SSS??^ - 1171?? + 352) \E\'^ 

+ (Xl+X2)l ^5 + p 

+ l{-27i' + —^ ^+538)\E\' + ^ ^ } 



2 J 

61 j 8r? (I3;r + 125;/ - .-xSl) |£l 4t/ (4;-!;r + 31!);/ - .192) jEp 2?] (lO;/- - 131;/ + 1008) 

+ 277(23277- 1097) I (27) 

3.3. The Elements of the Kepler equation 

_ 1 j v{v + 4)^\E\^-2L^\E\^ \ 
^~c4| 2V2L j 

Ss \ e^\E\^l\{13n - 8)(xi + X2) - (?7 + 8)vT^(xi - X2)) 1 
c4 I 2V2L^ J 

1 eN / XiX2(r;-4)7;|ii;|3/2 

+ ?^^Z3| 7i 

^^2|p|3/2/^^ A8 - 377)vT^ , (7/(27/ -19) + 8) ^^ 377VT^ , (7/^6)77 

. X^|.|- (C.2 (^^^^^ + i^fcM) . ^ . ^) ) 

1 f (I728e^ + 77 (33r/2 + 600r? + 37r2 - 4148)) L77 (237/2 - 477 - 64) \E^I'^ 

48\/2L3eN 4V2eN 

r] (105772 + 6277/ + 37r2 - 4232) |£;|5/2 1 

24y2LeN J 

^sj / v'2 - 877 (-19772 + 29287/ - 2976) \Efl'^ 

+ c6 jw-X2)(^ 64L4^N 

V2-87/ (-12772 - 145377 + 1544) \Efl'^ V2-8r/ (57772 + 8r/ - 256) |i;|^/2 
16Z/2eN 16eN 
/ (1307/3 + 23 597/2 - 72327/ + 2976) \Efl'^ (8677^ + 12027?2 - 37II77 + 1544) \Efl'^ 

+ (xi + X2) A^^. ^- 



(5O773 + 35772 - 47277 + 256) |£;|^/2 \ 1 / 377 (26772 + 10377 - 828) \Efl' 



8A/2eN )\^c^Y'^\ 8V2L^e^ 

7/(71772 + 36377 - 2388) |i^;|5/2 77(2772+977 + 52) |i^|V2 x | 

4y2X3eN 2^/2Le^ Jj' ^ ' 

^ 1 / e2 7/ (6772 + 1277 + 23) \E\^/^ 1 J| f xiX2eg,77 (27?2 - 357? + 4) \E\^/^ \ 
c6 \ 8x/2L3 J 4x/2i5 j 



/2 



2u 



5, f 4 V 5 - 2r? (I3r?^ + 327/ - 8) |E|3/2(xi - X2) 
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egj (2l7?3 + 23r?^ - 48r? + 8) \E\y^xi + X2) \ 

^ _ ^ , . , ^. e^vi2r,'-llr, + 4)\E\^/' 



48\/2L3 2AV2L J |' 

eN/? (27?2 _ llry + 4) |£:|3/2' 



8V2L^ 



+ -X2) 



+ — W5Z^ — ^ — W5Z^ — ) 



_ 1 372(5 -2r;)|j^|3/2 



+ |i5|3/2y2| '^5^ 4^/^^(r?-3)(xl-X2) 2 (r?^ - 8?? + 6) (xi + X2) | ^ ^| [ 2x1X2(15 ■ 
, Xi {Cq2 (Sy^ - '-f + 7)+VT^{CQ, (f - 7) + ^ - 4) + 2/?^ - + 4) 

, xf(gQi(3r/^^f^ + 7)+Vr^(CQi(7-f)-f^+4)+2^2-gfl+4) ]] 



1 / (l440r/2 + (I237r2 - 13952) r] + IOO8O) \E\^/^ 3 (ISt?^ - 55?? + 95) {El^/"^ 



\ 48\/2i3 2V2L 



Sg( /y2^(l3r?2-57r? + 75)|iJ|5/2 3 - 2r; (7r;2 - 99r; + 140) |i?|3/2 
+ ^|(Xi - X2) j^, 

, / 3 (2r;3 - 69tj^ + 239r/ - 140) \E\^/^ (-IS??^ + 1317^2 - 294r; + 150) \E\^^^ \ 
+ (Xi + X2) ^ + ■ ' ' 



5%i ( 3^27; (7??2 + 83?7- 315) |i;|3/2 r/ (77^ _ l98r? + 429) |£;|5/2\ 1 

+ [ T^Zi J I ' 

3.4.. The Elements of the Orbital phase 

1 \i.el{n-^^'')\^ 1 . {^h^seln{xi (??-x/r^-l)+X2 (tj + x/T^- l)) 1 
= ?\ MA / + ^^'\ AIA / 

6% { Xl (i^C-Qie^ {-Qrf - 19/? + 6) + {^e<^QieUG ' ^v) + ^) + is^U^ ' 2'?)^) 



+ 



X2 ' (i^CQ2e2, (-677^ - 19r? + 6) + ^/T^ (3^CQ2e^(77? - 6) - ^) + 3^,62,(3 - 27?)7?) 



L6 



_^ Xi X2 e^7y(277 - 5) 



4L6 j 
1 f -4O773 - 3847?2 + (497r2 - 1076) 7? + 256 
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_^ (-807?^ + 336r?2 + (1192 - 4977^) rj - 256) |£;| _^ rj {2Arf - 40r? - ll) \E\'^ 



128L4 16X2 



5, 



sj ^/^ ^r?(-212r?2 + 119,; + 492)-4e 
+ (Xi +X2)I ^ 



+ 



T] [1527]^ - 47r/ - 296) \E\ 

77 (I04r;^ - 1277? + 24) \E\^ \ _ / ^vT^^r;(31977 + 172) - Ae^y/T^ 

16X3 I + (Xl X2 ) I ^7 

_ VI - 47/r?(197?7 + 136)|£;| 3V1 - 47?j?(37r? - 8)\E\A 1 

leZs + 16Z3 J I 

si j fv (18()//- + 89.-.// + 47(i) ;/ (80;/- + :-!79// + 2(i()) \E\ 
r?(2?72 + 137r?-46) I^P 

1 3egj?7^ '^S / egj(7?- l)77(xi +X2) 4\/l - 4????(xi - X2 ) ( 1 .2 J , Xi X2 - 4)?? 



+ 



+ 



XI ^ (-^CQiegjr?(2r? + 1) + ^/T^ - ^Cqieg,/?) - ^6^(7? - 2)7?) 
X2 ^ (-^CQ2egj7?(27? + 1) + y/T^ (^gQ2egjr? - ^fe') - ^6^(7; - 2)77) 

Z6 



1 j 7^ (220 + 377^) eg,7? + a^eN7?^(157? + 32) _ eN(257/ + 52)7?^|.E| eN(2677 - 9)772|Sp 



c6 1 X6 64L4 64L2 

/ eN7? (-17077^ + 77; + 96) eN7/ (1887?^ + 437? - 144) 

"'"„fil(.Xl+X2jl loor? "I" 



.6 yyA.. ^28X7 64X5 

eN7?(l57/2 - 17?7 + 3) , /3eNVr^4^?7(6l77 + 32) 

+ (.Xl ~ X2J 



16X3 J -^^'V 128X7 

_ 3eNx/r^^??(597? + 48)|X;| _ 3eN(l - 4?7)3/27y|x;| 
64X5 16X3 

Si j le%r] (3772 - 77 - 6) + ^eN7? (677^ + 3377 - 2) 

+ ^XiX2| ^8 

eN7? (5077^ + 211?7 + 32) |^| eN?7 (57?^ - I777 + 4) \ 

32X6 32X4 J ' 

1/7/(10772 + 2877 + 5) 77 (IO772 + 2877 + 5) 77(107/2 + 28?7 + 5) |£;|' 



6 128X6 32X4 32^2 



2 / 77 (12772 - 1377 + 20) 77 (12772 - 1377 + 20) 77 (12772 - 137/ + 20) |x;|^ 

+ '^sXiX2 ^778 TTTTi^ + 



64X8 16^6 16X4 / 

7/(40772 + 3177-42) 77 (407/2 + 317/ - 42) |£;| 7/(40772 + 3177-42) |x;| 



2 



+ Ss< (Xl + X2) r— 77 + 



128X7 32X5 32X3 
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f 2lVl^r,{r, + 2) 21^1^7^(7^ + 2)\E\ 2lVl^7^{fj + 2)\E\^ \ 
+ (Xi-X2)(^ 128i^ 32L^ + 321? 1^^' ^^^^ 

I if3{xi-X2)v'VT^{l-2L^\E\f' 3vH2v- I) {X1+X2) {I ~2L-'\E\y 

^5v — -(^■{OS 



c6 1 \ 128i7 128i7 

,5/2 ^ 



5773(1 -2X21^1)- 



(35) 



As the binary system looses energy and angular momentum via the emission of gravitational waves, the above orbital 
elements will not remain constants of motion. One can evolve the binding energy and the angular momentum via 
balance equations connecting the far-zone flux with the near-zone, deduced from the time derivatives of the source 
multipole moments. This will be done in the following section. 



4. Energy and Angular Momentum Loss 

For the radiative dynamics, 2PN point-mass contributions and the leading-oder terms for the spin-orbit, spin(l)- 
spin(2), spin-squared interactions and tail contributions [34] are used. We work in ADM coordinates and use 
canonical spin variables (see e.g. [35] and references therein). 



4-.1. Source multipole moments in harmonic coordinates, using covariant SSC 



We begin this section by collecting the relevant source terms of the far-zone gravitational field. It is to be mentioned 
that we can rewrite spin contributions in terms of spin tensors rather than in terms of spin vectors. 



Sr, ^ S, 



a {k)(r 



The inversion is realised via 
The far-zone field of the gravitational wave reads 



= {1,2} (36) 



(37) 

defining TZ to be the distance observer-center of mass of the 



binary, and Na, NaiNa2---Nai to be an Z-fold tensor product of the line-of-sight vector Na^, 



lTT 




1 (3) 4 (2) 



(4) 



(3) 



(5) 



(4) 



60 



^abclijahc ~^ ^^^ah{iJ j)acd ^bcd 



^t^^ablijah c2^ab{iJ j)Q 



(6) 



be 



(5) 



^QQ-^abcdlijabcd ~^ ^Q^ab{iJj)acde ^bcde 



(38) 



The point-mass multipoles appearing above at 2PN are taken from [22], the linear-in-spin terms from [37], and 
the leading order spin-squared contribution from [24]. Note that [22] applied the harmonic gauge and [37] the 
covariant spin supplementary condition. In the following, is regarded as a bookkeeping parameter. To check out 
the relative orders of magnitude, we employ the scaling according to Sec. Appendix B and look for the powers of 
inverse c in such a way that the units of spinless an spin dependent contributions are the same. We remind of the 
fact that a time derivative goes along with a factor . Written in terms of spin tensors satisfying the covariant 
spin supplementary condition, they read 



1 

42 c2 



(29 - 8777)^2 -(30- 48r;)- 



1 



-^(253 - 183577 + 354577^)?;'* -t- -^(2021 - 594777 - 488377^)- v 
504 756 r 



2\ ^ ..2 



— (131 - 90777 + 127377^)- 7^2 
756 r 



1 
252 



(355 + 19067/ - 33777^)- 



42 c- 



■(24-7277) 



rr 
7a 



(26 - 202r/ + 4187;2)v2 
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sisi 



Jij 



Hjk 



Jijkl 
Jij kin 



+ -^(1085 - 405777 - 1463772)- 



— (11 - 3377) + ^' 



5 



63 
^s4 



21 c2 

(1 - 577 + 5r/2)r-2 



(41 - 33777 + 73377^)7;^ 



126 



^(742 - 3357, - 9857,2)1 



+ ^ -^ri'G-'n? STF,,- {^XiV^^ - v^x"^) (s^ + ^2 (,)(,)) 
+ ^-STFi.G^m^ f Cqi ((x/T^ + 3) r, - - l) ~S\S{ 

+cq2 (- (yr^ - 3) 77 + v/r^ - 1) 5^5^) I , 

S\ {m){k)S\ (m)(k) — S\ (i)(k)S\ {j){k) , 



r]^/^^G^m^STFijJ x'^^ [l + ^((5 - 197?)7;2 - (5 - 1377)^) 



- x'^v'' 



rr 



(1 - 277) 



+ x^'v^'' 



(1-377)770^5 STFijkAx 



r(l- 277) 



(1 - 377) 



+ 



110 C2 



(103 - 73577 + 139577^)7;^ - (100 - 6IO77 + 105077^) 



v'x^''^ 



72 rf 
55 c2 



(1 - 577 + 577^) \ +v''^x' 



78 

55 c2 



(1 - 577 + 577^) 



lijkim = ^^/^^77(277 - l)G^m^STFijkim {cc'^'"'"} , 



- 4:r]r]G^m^STFijejab\ x'^v^ fl + — ((13 - 6877)7;^ + (54 + 607,)^ 



+ v^x" 



(5,9 3 



rr 
28c2 



(5 - lOr,) 



+ ^ - 77G2 7n3 STF,,- { x'e^''^ + l)S^ (k)(m) + - 1)^2 



{k){m) 



= ^77 (37,- l)G3m4STFi,-;i [ekabX^^^'v'] + ^iiG^m*Si:F,,k\^-'^ {^rf - 577 + l) 



45 
2 

9r 



(577^ -577 + 1) r'^ekmnv'v^v^'x'' 
(43772 + 877-7) ek^nv"'x''x'x^ 



- ^ (92577^ - 38577 + 41) e,,mnx'v'^x^xh''^ , 
lyr^(l - 277) 77 GWSTF,,fe, {e/a6X»'^ V} , 
4i»7 (1 - 577 + 577^) 771^ STF,,fc,„ {e„a6 x'^^^'^'t;"} 
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where the symbol STF^^ . i^. labels the symmetric and trace-free part of an expression with respect to its indices 
ii,. . ■ ,ik- We do have access to higher-order accelerations in the ADMTT gauge, as we use Poisson brackets to 
compute time derivatives of separations and spins. We need a contact transformation that arbitrates between the 
covariant spin supplementary condition and canonical spin variables on the one hand and between harmonic and 
ADM gauge on the other. We give the relevant transformation terms in the subsequent section. 



4-.2. Transformation from harmonic to ADM and from covariant to canonical spin variables 



The coordinates transform from ADM to harmonic gauge (subscript "A") according to [22], and from the canonical 
spin variables to covariant spin supplementary condition (subscript "HC") as in, e.g. [37, 4], 

1 



'-HC 



2(t){k)) 



+ I {-V rl rl + 5r]rlvl 24 ta + 2rA) - 18 Tjrl v^rA ta } , 



1 

da; 



(48) 
(49) 



^HC 



"-HC 



dt 



HC 



dx 



HC 



dtf 



c2 2r2 



-»^12A 



"-12A V 



dtA 
d^HC 

3fi - 7v 

TA 



■>2{i)(k) 

38 



rA + v^lv 



9rl 



5vl 



7'A 



34 




(50) 



In the appendix of [1], the transformation of v was incorrectly implemented at 2PN point mass level. The gauge 
term in Eq. (A2) of [38] has to be reproduced, and finally reads 



g^(0+2)PP+PP 
Sx.+ 



V ■ 



I nn 
/ nv 

2r2 ' 



(t] (-55rf^ -I- Srw^ + 24) + 2) 



(^7(3 



T,(X, + ) 

Irv"^ - 20) - 4) + ^'"^^ (?7 (9rf2 



5rv^ + 34) + 2) 



(51) 



The operator V projects onto the transverse traceless part of the radiation field hij , see Section 7 of [1] , and for 
convenience, the superscript (0-|-2)PP+PP stands for the 2PN point particle contribution solely of the leading-order 
gravitational wave polarisation ^_|_.x with velocities and distances modified according to Eqs. (48) and (49). All 
the distances and velocities in Eqs. (39)-(47) have to be understood as harmonic variables with covariant spin 
supplementary condition, and we express the results for energy and angular momentum loss in ADM coordinates 
from now on for the remaining sections of the paper. 

The next section applies expressions of the far-zone flux of energy and angular momentum of the binary and 
provides differential equations of the binding energy 



4-3. Differential equations for the orbital elements 



In this section, we derive (orbital averaged) differential equations for the loss of orbital energy and angular 
momentum in terms of the energy itself and the radial eccentricity e^.. 

The balance equations for the decay of Li and \E\ through 2PN read [20] 



d\E\ 
dt 



G 



1 (3) (3) 



1 



1 (4) (4) 
lij k lij k 



189 



16 (3) (3) 

45 



1 



1 (5) 

9072 



(5) 

lijkrn I ijkin 



1 (4) (4) ■ 
^ JijkJ ijk 



(52) 



(dL, 
V dt 



G 



2 (2) (3) 1 
- Ipjiqj +^ 



1 (3) (4) 



32 (2) (3) 



63 
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where a superscript in round brackets, (n), denotes the n"^ time derivative. Note: the spin-squared contribution 
to are static in the sense that they contain neither velocities nor distances, and as one inserts the equations of 
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motion it becomes clear that they do not contribute. Without specifying the direction of the spins, the results for 
the instantaneous parts read 
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The tail terms at leading order can be taken from [34]. We use the quasi-Keplerian parameterisation from the 
previous sections to express the time dependent terms in terms of the eccentric anomaly £ when we specify to 
aligned spin vectors and orbital angular momentum. The time average of energy loss is done using the following 
relation, 

r£=2Tv 



d\E\ 



d\E\ 



dt / ^ Jg^o dt 



df 



df, 



(56) 



where the term in square brackets results from the Kepler equation. We can insert the above integration limits 
because the Kepler equation possesses fixed points at £ = mr, (n G Z). The squares of the velocity and the radial 
velocity as functions of £ enter via 

=r'^ + r^<i)^ . (58) 

The orbital averaged differential equations show total agreement up to all orders (purely including the 
instantaneous parts) appearing in [22], and symbolically read 
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(61) 



The symbol (. . .) means average over one orbital (radial) time period. The spin-orbit terms match with [39]; the 
spin(l)-spin(2) terms match [40], where we checked the aligned spins case. The explicit expressions read 
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5. Conclusions 



In this article, we have completed a previous work on (anti-)ahgned spins in a compact binary system [1] by next- 
to-next-to-leading order spin effects and by spin dependent radiation reaction effects. We provided expressions for 
the decay of a set of orbital elements, namely the binding energy and the radial eccentricity (|i?|,er); the reader 
may rewrite the above expressions to other sets of integrals, e.g. to use {Af,et) instead of {\E\,er). The results 
for spin-squared far-zone flux at leading order show up to be conform with [41], where the term of interest for the 
binding energy come purely from the equations of motion and the higher-order multipole moments rather that the 
leading-order quadrupole term in the mass quadrupole. The 3PN point-mass contributions to the energy loss are not 
included, but can be taken directly from the literature for a further publication, as well as the spin contributions, 
as soon as they are available. 

It is up to the reader to complete the discussion about conservation of alignment by proving conservation of 
constraints in [1] to higher orders in spin-orbit and spin-spin interactions, where in the case each spin vector points 
in the same direction as L nothing has to be shown. 

A subsequent publication will discuss the approximate solutions to the conservative dynamics in general orbits 
and arbitrarily orientated spin axes under the influence of spin-spin and spin-orbit interactions. For a naive insight, 
one could - for simplicity - assume that the time of observation of the binary is of the order of several orbital 
revolutions (rather than the much larger spin precession time scales) and assume the spin dependent terms to be 
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approximately constant. Then one is able to employ the methods for a calculation of eccentric orbits from the 
literature which have been used here. 
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Appendix A. Meaning of symbols 

The labeling is according to [42]. For convenience of the user, we list the complete set of used symbols in the present 
article §. 

quantity description 

c"^ power counting for post-Newtonian orders, mostly set to 1 

nPN n'^ post-Newtonian order, ©(c"^") 

Ss spin power counting, mostly set to 1 

X a quantity related to orbital angular velocity 

\E\ absolute value of binding energy 

\L\ angular momentum of orbit, L :— |L| 

Sa spin vector of particle a 

Gqa quadrupole constants 

e^x.y.z} orthonormal basis in local Lorentz frame 

Xa projection of spin a onto e^: Xa := (Sa • e^) 

e{j ,,,0} "time", "radial", and "phase" eccentricity 

ttr semimajor axis 

Tf symmetric mass ratio: 77 := 17111712 / {mi + m2)^ 

M. mean anomaly 

Af mean motion or radial angular velocity, respectively 

£ eccentric anomaly 

V true anomaly 

dynamical phase as function of £ 

$ orbital phase between two consecutive periastron passages 

Zj„ mass-type multipole moments 

current-type multipole moments 



Appendix B. Dimensionless Quantities 

Everything appearing in our prescription and the code is evaluated in scaled (dimensionless) quantities. The 
scaling is as follows: 



Gto- 



(B.l) 
(B.2) 



t 



E 



_G7n 



(B.3) 



r 




(B.4) 



P 



pc^. 



(B.5) 



§ Note that (in contrast to [1]) the orbital angular momentum is now denoted by L. 
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Si=Si^, (B.6) 

c 

S2 = S2^, (B.7) 
c 

where oj = d(j)/dt is the (unsealed) angular velocity. The bars are, from now on, omitted: "bared" quantities on 
the right-hand sides are understood as to be used in each case. \E\, which will be mostly used later as the binding 
energy will increase due to radiation reaction. 
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